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Abstract 

We use the non-perturbative QQ potential at finite temperatures derived in the 
Field Correlator Method to obtain binding energies for the lowest eigenstates in 
the QQ and QQQ systems (Q = c,b). The three-quark problem is solved by the 
, hyperspherical method. The solution provides an estimate of the melting tempera- 

ture and the radii for the different diquark and triquark bound states. In particular 
we find that J/tp and ccc ground states survive up to T ~ 1.3 T c , where T c is the 

^— > ■ 

critical temperature, while the corresponding bottomonium states survive even up 
to higher temperature, T ~ 2.3T C . 
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PACS numbers: 12.38Lg, 14.20Lq, 25.75Mq 



1. Introduction 



(N 

Diquark and triquark correlations in the quark-gluon plasma (QGP) are important for 
understanding the dynamics of heavy-ion collisions, processes in the early Universe and 
possibly the cores of the neutron stars [lj. Whether hadrons survive in the deconfined 
QGP is one of the key questions in QCD. In particular, the behavior of charmonia across 



O. — ' n 

the deconfinement transition has been suggested □ as hard probes of the QGP. Above T c , 

the critical temperature, most part of static QQ interactions disappear that eventually 
. J— i , 

^ . implies a dissolution of heavy quarkonia bound states into the continuum. A suppression 
of heavy quarkonia production in heavy ion collisions is usually considered as an observed 
signal of deconfinement. 

On quite general grounds it is expected that the QQ interactions get modified by tem- 
perature. As it was first mentioned in Ref. [3] an important part of the static interaction 
survives at the temperature above transition and can support the bound QQ states even 
after the deconfinement transition. Later existing of charmonia bound states at T > T c 
was confirmed on the lattice [I]. 



1 This effect was first first investigated in a phenomenological potential model taking into account the 
Debye screening [2]. 
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Most recently, in line with the study of Ref. [3], a new approach called the Field 
Correlator Method (FCM), was proposed to study dynamics of QGP, where the main 
emphasis was done on the nonperturbative vaccuum fields [S] , [E] , [7] . Let us summarize 
the basic formulation of the FCM as applied to finite T. The approach is based on the 
study of the quadratic field correlators < tr F^^&i^x, 0)F\ a (0) > (x is Euclidian), where 
0) is the parallel transporter necessary to maintain gauge invariance. The Gaussian 
correlator is expressed in terms of two scalar functions, D(x) and Di(x), which define the 
static potential between heavy quarks at T = 0, the confinement part of the potential, 
ar, of is expressed only in terms of D(x): 



a 



oo oo 

2 J d\ J drD(V\ 2 + r 2 ). (1) 



At T > T c one should distinguish between electric and magnetic correlators D E (x), 
D H (x), D E (x), and D±(x), and, correspondingly, between a E and a H . It was argued in 
[8] and later confirmed on the lattice [9] that above the deconfinement region D E (x) and, 
correspondingly, a E vanish, while the colorelectric correlator Df (x) and colormagnetic 
correlators D H (x) and D±(x) should stay unchanged at least up to T ~ 2T C . However, 
the correlators D H (x) and D^(x) do not produce static quark-antiquark potentials, they 
only define the spatial string tension a s = a H and the Debye mass vtld oc ^fo~ s that grows 
with the temperature in the dimensionally reduced limit [10J. 

The main source of the quark-antiquark static interaction at T > T c originates from 
the correlator D E {x) pj] 

l/T r 

V£*(r,T) = J du(l-uT) J \d\D Einp \VW + ~^), (2) 
o o 

that is responsible for bound states of quarks and gluons in the QGP . 

The purpose of this Letter is to report our latest results of low lying charmonium and 
bottomonium as well as ccc, bbb bound states above T c , obtained from the use of this 
approach. The paper is organized as follows. In section E] we introduce the static QQ 
potential. In section ED the numerical results for the ground states of J/ip and T are 
presented. In section ^} we evaluate the masses of the three-quark ccc and bbb baryons. 
Our conclusions are given in section [5] 
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2. The heavy quark-antiquark potential at T > T c 

In the framework of the FCM, the finite temperature behavior of the static QQ potential 
V Qg( r ' T) akT > T c is given by 

V QQ (r,T) = V^%T)+V^(r,T). (3) 

In Eq. §3$) Vq% (r, T) describes the interaction at short distances whereas V^(r, T) is the 
long-distance potential. The short distance interaction is represented by the perutrbative 
one-gluon exchange potential that undergoes a Debye screening by the color charges of 
the QGP 

V^\r, T) = -C F ^ exp(-m D r), (4) 

where Cp = 4/3 is the color factor, is the inverse of the Debye screening radius. This 
form of the modified Coulomb potential has been used in earlier works (see e.g. Ref. [T2] ) 
to specify the in-medium potential between heavy quarks and determine the dissociation 
points of different quarkonium states. The Debye mass in Eq. (j4J) can be written 
as rriD = 2.06 v/cr s (T), where <J S (T) is the spatial string tension due to chromomagnetic 
confinement. In what follows we use the results of Ref. p2] where the quantity yJa s (T) 
has been computed nonperturbatively up to two loops in the deconfined phase of QCD. 
For rif = the Debye mass varies between 0.8 GeV and 1.4 GeV, when T varies between 
T c = 275 MeV and 2 T c . 

The long-distance interaction V np (r,T) requires theoretical assumptions about its 
shape. We follow Ref. [Ej, where the long distance nonperturbative QQ potential was de- 
rived analytically from the analysis of the non-perturbative part of the correlator function 
Df(x): 

D? {np \x) = 5exp(-M |x|)/|x| (5) 

In Eq. (JHJ) the coefficient B must be considered as functions of the physical temperature. 
In the confinement region 

B = 2C F a f s a adj M = Qa f s a f M 0} (6) 

a{ being the freezing value of the strong coupling constant in the confinement region, <jf 
is zero temperature string tension, and Mq has the meaning of the lowest gluelump mass 
PS]. Taking a{ = 0.6, a f = 0.18 GeV 2 , M = 0.9 GeV, we get B = 0.583 GeV 3 . Above 
the critical temperature one substitutes B — > £(T)B, where the T-dependent constant 

e(T)= (l-0.36§^) (7) 
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Figure 1: The dependence of the V^(oo, T) on the temperature T given by Eq. (jHJ). 



arises from the comparison with lattice data [llj. As the result one obtains 

V^(r,T) = £(T) (y^(oo,T) - V^(r)) , 
where the continuum threshold (a constant shift in the potential) V np (oo,T) in Eq. 



is 



V np (oo,T) 



B 

Ml 



1 - 



T 



exp( 



Mo- 
T ■ 



and 



V np (r) = xK l {x) + 0[ — ). 



T 

Mo 



(9) 



(10) 



In Eq. ffTUj) Ki(x) is the modified Bessel function and x = M r. At T = T c one obtains 
V^(oo,T c ) ~ 0.5 GeV that agrees with estimate obtained in Ref. [H] from lattice data. 
The large positive value of the continuum threshold is a consequence of non perturbative 
vaccuum fields and can not be explained in perturbation theory. The behavior of the 
potential Vqq(oo,T) is shown in Fig. [T] 



3. Quarkonia at finite temperature 

Having specified the static QQ potential, we can now exploit the relativistic Hamiltonian 
technique, developed in [16] to calculate the masses of the QQ states as a function of the 
temperature. This technique does not take into account chiral degrees of freedom and 
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is applicable when spin-dependent interaction can be treated as perturbation. Therefore 
below consider heavy quarkonia and baryons. The masses of heavy quarkonia are defined 

as 

Win 

Mqq = + + £ K?,^q), (11) 

where E (rriQ, /xq) is an eigenvalue of the Hamiltonian H = H + Vqq, where Vqq is 

given by Eq. ([3]), trig are the bare quark masses, and [Iq are the constant auxiliary fields 

(AF) that were introduced to treat the kinematics of the relativistic particles. These 

parameters have to be found from the variational condition 

dM n n , , 

(12) 



9^Q 

The eigenvalue problem is solved for each set of hq; then one has to minimize Mqq 
with respect to fiQ. Such an approach allows for a very transparent interpretation of AF: 
starting from bare quark masses tuq, we naturally arrive at the dynamical masses \iq that 
appear due to the interaction. The AF are treated as c-number variational parameters. 
The bound state exists if E (mQ, //q) < V np (oo,T). 

To find out whether the non perturbative interaction can support the bound states at 
T ~ T c we use the Bargmann condition 

oo 

n<I = » Q f \V np (r) + V^(r,T)\rdr = ^ (J^- + , (13) 
o 

where n is the number of the S'-wave bound states. Since ~ 1 GeV, B/Mq ~ 0.5 
GeV, and M ~ 1 GeV, we conclude that to support at least one bound state one needs 
Hq > 1 GeV, i.e. there is no bound states of light quarks, but cc and especially bb binding 
is possible. Moreover, the bottomonium spectrum should display a much larger number 
of bound states above T c . 

The solutions for the binding energy for the IS" states of charmonium and bottomonium 
are shown in Tables [H, EJ In these Tables we present the constituent quark masses \iq for 
cc and 66 El, the differences Eq = E — Vqq(oo), the mean squared radii \/r^~, and the 
masses Mqq for the cc and bb mesons. We employ m c — 1.4 GeV, m& = 4.8 GeV, and 
a s = 0.35. As expected, we obtain the weakly bound state at T = T c that disappears at 
T ~ 1.3 T c . The charmonium masses lie in the interval 3.2 - 3.3 GeV, that agrees with 
the results of Ref. [TT]. Note that immediately above T c the mass of the cc state is about 
0.2 GeV higher than that of J/ij). 



2 These masses are computed solely in terms of the bare quark masses m c and mi,, respectively. Note 
that, as in the confinement region [17] . the constituent masses [iq only slightly exceed bare quark masses 
tyiq that reflect smallness of the kinetic energies of heavy quarks. 
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Table 1: The IS* J/ip state above the deconfinement region. Mo = 0.9 GeV, B = 0.583 
GeV 3 , m c = 1.4 GeV, a s = 0.35. Masses and energies are given in units of GeV, the 
mean squared radius V< r 2 > in units of GeV -1 . 



T/T c 






V< r 2 > 


M c - C 


1 


1.451 


-0.019 


7.53 


3.291 


1.3 


1.419 


+ 0.006 


10.50 


3.186 



Table 2: The IS* T state above the deconfinement region. The notations are the same as 
in Table GQ m 6 = 4.8 GeV 



T/T c 




- V Q q(oo) 


V< r 2 > 




1 


4.984 


-0.300 


1.27 


9.815 


1.3 


4.950 


-0.183 


1.55 


9.802 


1.6 


4.915 


-0.095 


2.06 


9.783 


2.0 


4.863 


-0.021 


4.25 


9.742 


2.2 


4.832 


-0.003 


7.62 


9.712 


2.3 


4.818 


+ 0.001 


9.52 


9.694 
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As expected, the bb bound states remain intact up to the larger temperatures, 
T ~ 2.3 T c , see Table [21 The masses of the L = bottomonium lie in the interval 
9.7-9.8 GeV, about 0.2-0.3 GeV higher than 9.460 GeV, the mass of T(1S). 

The results for cc and bb bound states found from the FCM can be compared to the 
calculations based on phenomenological QQ potentials identified with the free energy 
measured on the lattice [IB], [HJ]- Our results for lS(J/ip) are qualitatively agree with 
those of Refs. [18] . [TP"] while our melting temperature for lS'(Y) is much smaller than 
T ~ (4 - 6) T c found in Ref. [19]. 



4. QQQ baryons at T > T c 

The three quark potential is given by 

Vqqq = lY,V QQ ( rij ,T), (14) 

i<j 

where | is the color factor. We solve the three quark Schrodinger equation using the 
hyperspherical method. The wave function of a QQQ baryon depends on the three-body 
Jacobi coordinates 




(i,j,k cyclic). There are three equivalent ways of introducing the Jacobi coordinates, 
which are related to each other by linear transformations with the Jacobian equal to 
unity. In what follows we omit the indices i and j. 

In terms of the Jacobi coordinates the three-quark kinetic energy operator H is written 

as 

i / d 2 d 2 \ i / d 2 5 d L 2 (n)\ 

where R is the six-dimensional hyperradius that is invariant under quark permutations, 
R 2 = P 2 + A 2 = f (r? 2 + r 2 2 3 + r 2 3l ) , (17) 
p = Rsin9, \ = Rcos9, 0<9<n/2, (18) 

^2 denotes five residuary angular coordinates, and L 2 (VL) is an angular operator 

o9 2 o9 sin^ cos 2 ^ 

whose eigenfunctions (the hyperspherical harmonics) satisfy 

L 2 (Q)Y [K] (9,n p ,n x ) = -K{K + A)Y [K] (9, n p , n A ), (20) 
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with K being the grand orbital momentum. 

The wave function ip(p,X) is written in a symbolical shorthand as 

[K] 



(21) 



where the set [K] is defined by the orbital momentum of the state and the symmetry 
properties. We truncate this set using the hypercentral approximation K = K min = 0. 
Introducing the reduced function u(R) 

1 u(R,T) 



*(p,\r) 



(22) 



<7r"3 #5/2 

and averaging the three-quark potential (tHJ) over the six- dimensional sphere one obtains 
the one-dimensional Schrodinger equation for u(R,T) 



d 2 u(R,T) 
dR 2 



+ 2 



E - - ^£(T) (V pert (i?,T) + V np (R,T)) 



where 



vpcrt/D rp\ 4 ac{R) 
V p {R,T) = --a a - . 



«(i2,T) = 0,(23) 



(24) 



a c (5) 



16 //XQ 



tt/2 



7T 



cxp 



-mo R sin 8 



sin cos 



(25) 



and 



V np (i?,T) = y np (oo,T) - 



165 

7tM V /Xq 



tt/2 



#1 



M 0J Rsin6» 



sin 3 cos 2 0d0 5.(26) 



The temperature dependent mass of the colorless QQQ states is defined as 

3 m o 3 
M QQQ = o + o^Q + Eo(m Q ,n Q ). 

The bound QQQ state exists if E (rriQ, hq) < Vqqq(oq,T), where 

V QQQ (oo,T) = ^ p (oc,T). 



(27) 



(28) 



In Tables |3] and H] we show the masses of the ground ccc and bbb states as a function of 
the temperature. 
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Table 3: The ground ccc state as a function of the temperature above the deconfmement 
region. The notations are the same as in Table Q] 



T/T c 




E ~ V QSQ (oo) 


V< R 2 > 


M ccc 


1 


1.474 


-0.046 


6.60 


4.922 


1.3 


1.434 


-0.002 


9.96 


4.769 



Table 4: The ground bbb state as a function of the temperature above the deconfmement 
region. The notations are the same as in Table [T] 



T/T c 




E - Vqqq(oo) 


V< R 2 > 


M bbb 


1 


4.977 


-0.506 


3.12 


14.665 


1.3 


4.927 


-0.314 


3.86 


14.661 


1.6 


4.885 


-0.171 


5.04 


14.645 


2.0 


4.840 


-0.046 


8.21 


14.598 


2.2 


4.822 


-0.011 


10.59 


14.562 


2.3 


4.815 


+ 0.001 


11.62 


14.541 
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5. Conclusions 



Based on the non-perturbative dynamics driven by the field correlators at finite tempera- 
tures we have calculated binding energies for the lowest eigenstates in the cc, ccc, 66,and 
bbb systems. For what concerns the charmonium we find that J/ip (and ccc ground states) 
survive up to T ~ 1.3 T c . On the other hand, the bb and bbb states survive even up to 
higher temperature, T ~ 2.3 T c . This suggests that the systems are strongly interacting 
above T c . 

This work was supported in part by RFBR Grants # 08-02-00657, # 08-02-00677, 
# 09-02-00629 and by the grant for scientific schools # NSh.4961.2008.2. 
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